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' Formaldehyde  decay  was  followed  by  IR  emissionUt  3.5  >inr  behind  reflected 
shocks  at  total  concentrations  near  2.4  and  4.5  x /fpv^molecule/ertl^  Data  were 
obtained  at  O.S-tfS  Intervals  over  a temperature  range  of  1700-2710  K for  mixtures 
containing  0.1%,  0.5%,  and  1.0%  CH2O  In  Ar.  Exponential  decay  was  observed  In 
all  cases;  at  lower  temperatures  there  was  evidence  of  an  Induction  time  prior 
to  the  onset  of  the  decay.  An  apparent  rate  was  defined  as  1^^*  -[d(1n  signal)/ 
dt]/[M];  ^was  observed  to  be  first  order  In  argon  and  to  decrease  with  Increas- 
ing dilution.  yvt  t 

These  data  were  analyzed  In  terms  of  the  mechanism 

VCHgO  + M = HCO  + H + (1); 

^ " 2 U);  . 



VhCO  + H * H2  + CO^  (4) . 

Using  literature  values  for  k2,  k^,  and  k^,  a reasonable  fit  to  all  the 
data,  both  In  terms  of  temperature  and  concentration  dependence,  was  obtained 
with  l(^)=  6.0  x(l0^exp(-364  kJ/RT)  ^i^'^nolecule^ s^.  Evidence  Is  also 

presented  to  suggest  that  the  reactiof^^ Tki  ~ 

C:>  CH2O  + M = CO  + H2  + 

does  not  play  a significant  role  In  the  high  temperature  pyrolysis  of  formaldehyde. 


' CH2O  + H = HCO  + H2  ^ 
'ShCO  tM-H  + CO  + H 
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INTRODUCTION 


Formaldehyde  has  long  been  recognized  as  an  Important  intermediate  in  many 

oxidation  systems.  Furthermore,  several  recent  studies  have  suggested  that  the 

dissociation  of  formaldehyde  is  a particularly  important  reaction  in  such  systems. 

For  example,  Peeters  and  Mahnen^  reported  this  reaction  to  be  rate-determining 

in  CH^/02  flames.  In  a similar  vein,  the  sensitivity  analysis  of  Boni  and  Pen- 
2 

ner  indicated  that  formaldehyde  dissociation  was  one  of  the  four  most  important 
reactions  in  a typical  shock  tube  study  of  methane  oxidation.  In  spite  of  its 
importance,  there  have  been  only  two  high  temperature  studies  which  have  focused 
directly  upon  this  reaction.  The  first,  by  Gay,  et  al.,  measured  formaldehyde 
decay  in  2-5%  mixtures  with  neon  as  diluent  between  1530-2100  K.  In  these  ex- 
periments exponential  decay  of  formaldehyde  was  observed  after  an  induction  period. 
The  reaction  appeared  to  be  first  order  in  neon,  and  the  activation  energy  for  the 
observed  rate  constant  was  116±10  kj/mole.  Isotopic  scrambling  was  observed  when 
equimolar  mixtures  of  CH2O  and  CD2O  were  used.  This  observation,  coupled  with  the 
presence  of  an  induction  period,  led  these  workers  to  suggest  a variety  of  mechan- 
isms involving  hydrogen  atom  intermediates.  One  such  mechanism  was: 


CH2O  + M = H + HCO  + M (1 ) 
CH2O  + H = H2  + HCO  (2) 
HCO  + M = H + CO  + M (3) 
HCO  + H = H2  + CO  (4) 


The  main  difficulty  here  was  that  this  mechanism,  analyzed  via  steady-state  as- 
sumptions, predicted  an  activation  energy  for  the  observed  rate  constant  in  ex- 
cess of  209  kJ/mole,  much  higher  than  that  observed.  Other  mechanisms  proposed 
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also  had  certain  objectionable  features,  and  It  was  not  possible  to  assign  rate 
constants  to  the  various  elementary  reactions.  Probably  the  most  salient  feature 
of  this  investigation  was  the  convincing  evidence  that  secondary  reactions  played 
a major  role.  Thus  the  rate  constant  for  reaction  (1)  was  appreciably  smaller 
than  the  observed  rate  constant  for  CHgO  decay. 

4 

Schecker  and  Jost  studied  the  pyrolysis  in  argon  diluent  between  1400-2200  K. 

Here  the  distinctive  feature  was  use  of  a variety  of  CHgO  mole  fractions;  these 

ranged  from  0.02%  up  to  1%.  Qualitatively,  these  data  were  similar  to  the  earlier 

work— exponential  decay  with  induction  periods  at  lower  temperatures.  The  active 

tion  energy  was  seen  to  increase  from  144  kJ/mole  to  230  kJ/mole  as  the  CHgO 

mole  fraction  decreased  from  1%  to  0.02%.  It  is  important  to  note  these  two 

shock  tube  studies,  with  completely  different  diagnostics,  yielded  raw  data  in 

reasonably  good  agreement  with  each  other.  For  example,  at  1800  K,  Gay,  et  al., 

reported  an  apparent  rate  constant  of  3.5  x 10’^^  cm^  molecule”^  s"^  for  their 

-15 " 

2-5%  mixtures,  whereas  Schecker  and  Jost  reported  a value  of  1.2  x 10  for 
their  1%  mixture.  The  dilution  study  would  suggest  that  a somewhat  higher  value 
would  be  expected  for  the  more  concentrated  case.  Furthermore,  the  Schecker- 
Jost  values  were  obtained  at  ten  to  eighty  times  higher  total  concentrations. 

Thus  part  of  the  remaining  discrepancy  may  be  due  to  the  fact  that,  at  the  higher 
concentrations,  reaction  (1)  is  no  longer  in  the  low-pressure  region  for  a uni- 
molecular  dissociation  rtiction.  If  this  were  the  case,  the  primary  dissociation 
would  not  be  quite  as  rapid  as  calculated  at  the  higher  pressures;  this  could  be 
reflected  in  a lower  apparent  decay  rate.  In  any  event,  it  would  appear  the 


data  is  reasonably  consistent. 

Schecker  and  Jost  analyzed  their  data  in  terms  of  the  mechanism  given  above 


and  reported  values  for  the  various  rate  constants.  The  one  troubling  aspect 
of  this  analysis  Is  the  observation  that  the  activation  energy  for  the  1%  mix- 
tures, even  though  30  kJ/mo1e  higher  than  that  reported  by  the  earlier  workers. 

Is  still  quite  low  to  be  consistent  with  the  results  of  a steady  state  analysis 

* 

of  the  mechanism  listed  above. 

The  situation  with  respect  to  formaldehyde  dissociation  became  considerably 
more  confused  as  a result  of  the  CH^/Og  flame  studies  of  Peeters  and  Mahnen.^ 

Here  the  data  were  interpreted  to  indicate  that  formaldehyde  dissociated  directly 
to  molecular  products: 

CHgO  + M -V  CO  + Hg  + M (!') 

with  k^'  = 3.5  X 10'®  exp(-146  kJ/RT)  cm®  molecule"^  s'^  . In  addition  to  the 
obvious  difference  in  the  mechanism  of  dissociation,  this  study  suggests  a dis- 
sociation rate  at  1800  K some  three  orders  of  magnitude  larger  than  the  total 
decay  rate  measured  in  the  shock  tube  work. 

The  present  study  was  initiated  to  attempt  to  clarify  the  formaldehyde  de- 
cay process.  Here  an  extensive  numerical  integration  study  was  coupled  to  the 
traditional  data  reduction  techniques.  Thus  the  effects  of  the  secondary  reac- 
tions could  be  determined  without  resorting  to  approximate  schemes  such  as  the 
steady-state  assumption.  Furthermore,  this  approach  allows  one  to  properly  account 
for  the  temperature  changes  that  invariably  occur  in  such  a system.  It  was 
also  possible  to  extend  the  temperature  upwards  to  2700  K.  Analysis  of  this  new 
data  in  this  way  suggests  that  the  scheme  outlined  by  Gay,  et  al . , is  quantita- 
tively consistent  with  the  observations,  but  that  the  rate  constant  assignments 
of  Schecker  and  Jost  must  be  modified.  Perhaps,  most  importantly,  this  study 
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presents  additional  evidence  that  reaction  (1  ) must  have  a rate  constant 
many  orders  of  magnitude  lower  than  that  reported  In  the  literature. 


EXPERIMENTAL 


The  7.6  cm  shock  tube,  gas  handling  system,  and  optical  configuration  have 

5 

been  described  previously.  Infrared  emissions  were  first  collimated  by  two 

« 

silts  1.5  mm  wide  and  5.0  mm  high  placed  50  mm  apart.  Emissions  then  passed 
through  either  a 3.50  ym  Interference  filter  (FWHM  = 0.12  y)  or  a filter  wheel 
set  at  3.50  ym  (FWHM  = 0.15  y)  and  were  focused  upon  a llquid-nitrogen-cooled 
Indium  antimonide  detector.  The  3.50  ym  emission  should  be  due  exclusively  to 
CHgO.  The  formyl  radical,  a proposed  Intermediate,  has  its  band  centered  near 
4.02  ym.^  The  detector-preamp  combination  had  a nominal  rise  time  of  1 ys. 

Output  from  the  preamp  was  collected  by  a Blomatlon  Model  805  transient  recorder 
at  0.2  or  0.5  ys  intervals.  The  Interference  filter  was  used  for  the  most  dilute 
mixture;  the  larger  transmittance  here  meant  that  all  of  the  data  exhibited  simi- 
lar S/N  characteristics. 

Formaldehyde  was  prepared  by  heating  paraformaldehyde  to  'v-  100®  C,  passing 
the  distillate  through  a dry  ice-acetone  trap,  and  then  condensing  the  monomer 
at  77  K.  Care  was  taken  to  keep  the  monomer  pressure  below  1 kPa  (7.5  torr)  at 
all  times  to  minimize  polymerization.  The  monomer  was  then  subjected  to  several 
bulb-to-bulb  distillations  from  dry  ice  to  liquid  nitrogen  temperatures.  Then, 
while  at  -78®  C,  It  was  distilled  into  the  mixing  bulb.  In  most  cases  a fresh 
batch  of  CH2O  was  prepared  for  each  mixture.  The  argon  was  Airco  Research  Grade 
(99.9998%)  and  was  used  without  further  purification. 


The  shock  tube  was  typically  pumped  down  to  'v,  3 mPa  and  the  observed  leak- 
outgassing  rate  was  usually  near  4 mPa/min.  The  tube  was  Isolated  from  the  pumps 
for  approximately  one  minute  prior  to  shock  initiation;  the  nominal  background 
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pressure  was  near  7 mPa.  Since  the  test  section  was  pressurized  to  1.4-2. 7 
kPa  with  the  mixtures  of  interest,  the  background  impurities  from  the  shock 
tube  were  present  in  the  low  ppm  range  in  the  shocked  mixtures.  Mylar  dia- 
phragms and  helium  driver  gas  were  used  throughout. 

* 

Temperatures  behind  the  reflected  shock  were  computed  in  the  usual  way 
from  measured  incident  shock  velocities.  Reflected  shock  pressures  were  mea- 
sured with  a fast  response  pressure  transducer  and  were  always  found  to  be  in 
good  agreement  with  those  calculated.  Furthermore,  the  reflected  shock  pressure 
was  always  constant  over  the  time  interval  that  data  were  collected;  as  a result, 
no  attempt  was  made  to  correct  for  non-ideal  effects. 


RESULTS 

The  mixtures  studied  are  listed  in  Table  I.  Mixtures  A and  B were  prepared 

from  the  same  formaldehyde  batch,  whereas  all  the  other  mixtures  were  prepared 

* 

from  separate  batches.  In  most  cases,  the  3.50  ym  signal  decayed  exponentially 
after  an  initial  rapid  rise  upon  shock  front  passage.  At  lower  temperatures 
there  was  evidence  of  an  induction  period  prior  to  onset  of  decay,  and  this  phenome- 
non was  more  evident  for  the  most  dilute  mixture,  A typical  experimental  record 
is  shown  in  Fig.  1.  Each  experiment  was  characterized  in  terms  of  its  exponen- 
tial decay  by  an  apparent  rate  constant  k,  (k,  s - {d(ln  signal )/dt}/H)  where 

da 

M was  the  total  concentration.  Values  of  these  apparent  rate  constants  are  listed 
in  Table  I.  Note  in  the  Table  that  half  of  the  experiments  for  mixture  A were 
I done  at  total  concentrations  nearly  a factor  of  two  lower  than  all  the  other  shocks. 

I 

1 

Values  of  k^  there  fell  in  nicely  with  all  the  other  1%  data.  This  suggests  the 

decay  is  first  order  in  argon,  and  this  observation  is  consistent  with  that  of 

Gay,  et  al . , who  worked  at  comparable  number  densities.  All  of  the  data  with 
18  -3 

; M 4.5  X 10  molecules  cm  are  plotted  in  Fig.  2.  There  is  no  evidence  (cf. 

i 

Table  I)  of  variation  in  k^  with  different  mixtures  of  the  same  CH20/Ar  ratio, 
and  in  Fig.  2 mixtures  are  grouped  together.  Note  in  Fig.  2 the  decrease  in  k^ 

; as  the  system  becomes  more  dilute.  The  error  bars  shovm  are  for  one  standard 

deviation, and  it  is  evident  the  shift  is  not  simply  a statistical  aberration. 

This  trend  is  reinforced  by  another  sequence  of  experiments  with  a 0.2%  CH20/Ar 
mixture.  These  data  are  not  reported  due  to  larger  than  normal  temperature  un- 
certainties resulting  from  a malfunction  of  a velocity  gauge.  However,  these  data 
are  consistent  with  the  trend  seen,  particularly  with  respect  to  the  decrease  in 


L 
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k at  lower  temperatures  in  more  dilute  systems.  The  Arrhenius  parameters  cor- 

a 

responding  to  the  least-square  lines  in  Fig.  2 are  listed  in  Table  II.  Note  that 

the  trend  of  increasing  E,  upon  dilution  is  qualitatively  consistent  with  that 

observed  by  Schecker  and  Jost.  The  values  observed  here  for  k^  compare  favor- 

* 

ably  to  those  reported  by  the  other  shock  tube  workers;  at  1800  K,  the  present 
values  fall  between  those  reported  earlier.  The  most  significant  difference  here 
is  the  higher  value  of  E,  for  1%  mixtures;  note  this  value  is  over  50  kJ/mole 

a 

larger  than  that  reported  by  Schecker  and  Jost. 

i 


A 
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DISCUSSION 

Although  there  are  minor  differences  among  the  three  shock  tube  studies, 
as  a group  they  present  convincing  evidence  that  the  pyrolysis  of  formaldehyde 
is  a complex  process.  The  observed  induction  periods,  dilution  effects,  and  iso- 
topic scrambling  all  suggest  that  secondary  reactions  play  a major  role  here. 

A mechanism  consisting  of  reactions  (1)  - (4)  was  suggested  by  earlier  workers, 
and  various  attempts  were  made  to  check  the  consistency  of  this  scheme  with  the 
data  collected  in  this  investigation.  Following  the  earlier  workers,  application 
of  the  steady-state  assumption  to  both  hydrogen  atoms  and  formyl  radicals  yields: 

ka  = k^n  + (k2k3/k^k4)**]  [A] 

This  analysis  would  suggest  that  the  observed  activation  energy  should  lie  within 
the  range  (^)(E^  + ) (for  k2k2»k^k^)  and  (for  k2k2«k^k^).  For  the  for- 

mer  case,  it  was  assumed  that  the  activation  energies  of  k2  and  k^  were  negligibly 
small  relative  to  k^  and  k^.  Assuming  the  activation  energies  are  comparable  to 
the  bond  energies,  this  analysis  would  lead  to  an  observed  activation  energy  be- 
tween 'V  210  kJ/mole  and  'v  370  kJ/mole.  Although  the  experimental  values  are  com- 
patible with  this  lower  limit,  especially  when  it  is  recognized  that  k^  and  k^ 
may  have  somewhat  lower  activation  energies  since  these  dissociation  reactions 
are  occuring  in  the  low  pressure  region,  note  that  Eq.  [A]  predicts  that  there 
will  be  no  change  in  k upon  dilution.  This  prediction  is  in  marked  contrast 

a 

to  the  observations.  This  inconsistency  does  not  necessarily  invalidate  the 
mechanism;  the  problem  may  well  be  the  steady-state  approximations  used  in  the 
analysis. 

It  was  felt  a much  better  way  to  test  this  mechanism  would  be  to  numerically 
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Integrate  the  system  of  rate  equations  and  compare  the  calculated  profiles  to 

( 

those  observed.  In  addition  to  allowing  one  to  avoid  such  assumptions  as  the  j 

steady-state  approximation,  this  approach  allows  one  to  accurately  treat  the  ; 

I 

variation  In  temperature  which  accompanies  the  pyrolysis.  Neglect  of  such  varl- 

i 

atlons  could  be  particularly  troublesome  here  where  reaction  (1)  has  a large 
enthalpy  change;  the  problem  Is  exacerbated  by  the  varying  dilutions  used. 

The  procedure  employed  was  straightforward;  calculations,  using  a program  | 

previously  described,^  were  done  at  several  temperatures  for  the  various  mixtures.  I 

I 

! 

The  calculated  [CH20]-t1me  profiles  were  treated  In  the  same  manner  as  the  ex-  | 

perlmental  data.  Values  of  k.  calculated  for  a specified  Initial  temperature  could  • 

then  be  compared  to  those  observed.  Results  of  several  representative  calcula- 
tions are  shown  In  Table  III.  The  line  labeled  Cl  used  the  values  of  and  k^ 

4 8 

reported  by  Schecker  and  dost  ; k2  was  taken  from  the  recent  work  of  Kleitm, 

9 

and  k^  was  the  value  used  by  Bowman.  Note  that  the  calculated  k^  values  are 

much  too  high  here.  Results  of  attempts  to  lower  k^  by  using  the  much  lower  k^  ' 

g 

value  of  Bowman  are  shown  In  line  C2.  Note  such  a low  value  of  k^  completely  , 

suppresses  the  dilution  effect  seen  both  In  this  work  and  In  that  of  Schecker 

and  dost.  It  Is  evident  that  simple  variation  of  k«  Is  not  sufficient  to  recon-  ' 

die  the  calculated  and  observed  results.  , 

Of  the  three  remaining  rate  constants,  k,  would  appear  to  be  the  most  poorly 

• » 
8 

defined.  Reaction  (2)  has  been  carefully  studied,  at  least  at  lower  temperatures, 

and  reaction  (4)  Is  known  to  be  very  fast  even  at  room  temperature.^®  Equation 

[A]  suggests  k^  Is  at  least  as  sensitive  to  k-j  as  any  other  rate  constant.  In  ' 

this  light,  a series  of  calculations  were  done  where  k-j  was  varied  while  the 

other  rate  constants  were  fixed.  Here  the  Schecker-dost  value  was  used  for  kj 
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to  preserve  the  observed  dilution  effect.  This  value  is  similar  to  that  recently 

n 1? 

reported  by  Jachimowski , but  larger  than  that  obtained  by  Tsuboi.  A rea- 
sonable fit  to  all  the  data,  both  in  terms  of  temperature  and  concentration 
dependence,  was  obtained  with  = 6.0  x 10  exp(-364  kJ/RT)  cm  molecule  s“  . 

These  results  are  shown  in  line  C3  of  Table  III.  This  activation  energy  seems 
slightly  high  for  a unimolecular  dissociation  in  the  second  order  region.  Using 
a recent  value  for  the  formyl  radical,®  AEo  = 371  ± 10  kJ/mole  for  reaction 
(1).  Troe's  analysis  of  low  pressure  dissociation  reactions  would  suggest  an 
activation  energy  near  320  kJ/mole  for  reaction  (1)  with  this  AEo.  However,  use 
of  lower  activation  energies  for  k^  in  the  calculations  failed  to  give  agreement 
with  the  0.1%  mixture  data.  It  was  possible  to  use  lower  activation  energies  by 
increasing  the  temperature  dependence  of  the  secondary  reactions.  Here  the  most 
plausible  increase  involves  \f.2*  k^  already  has  an  activation  energy  equal 

to  its  endothermic ity  and  k^  is  expected  to  have  very  little  temperature  depen- 
dence. An  example  of  the  type  fit  that  could  be  obtained  is  shown  in  line  C4  of 
Table  III.  Note  that  a better  fit  could  have  been  obtained  with  a slightly  larger 
pre-exponential  term  for  k^,  but  such  fine  tuning  was  considered  unnecessary. 

This  adjusted  kg  is  equal  to  Klemm's  value  at  1800  K and  is  1.5  times  larger 
than  Klemm  at  2500  K.  It  appeared  that  no  real  improvement  was  afforded  by 
modifying  kg.  The  simpler  method  of  varying  only  k^  seemed  preferable,  parti- 
cularly in  light  of  the  obvious  uncertainty  in  k^. 

The  effects  of  varying  k^  and  k^  were  checked  by  separately  changing  these 
rate  constants  by  a factor  of  two.  Such  a change  caused  a shift  of  25%  in  k 

d 

for  the  1%  mixtures  at  1800  K.  For  other  conditions,  the  shift  was  smaller; 
particularly  significant  was  the  observation  that  these  changes  had  virtually 
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no  effect  upon  k.  for  the  0.1%  mixture.  Thus,  k,  only  depends  upon  the  values 

a a 

of  k-j  and  k2  here.  In  essence  then,  the  fitting  procedure  consisted  of  adjust- 
ing k^  to  fit  the  observations  on  the  0.1%  mixtures.  With  this  fit,  it  was  then 
found  that  good  agreement  could  be  obtained  for  the  more  concentrated  mixtures 
using  literature  values  of  the  other  rate  constants.  Furthermore,  given  the 
fact  that  k.  for  these  more  concentrated  mixtures  was  reasonably  sensitive  to 

a 

the  values  used  for  the  rate  constants  of  these  secondary  reactions,  the  fit 
obtained  suggests  these  literature  assignments  are  reasonably  accurate.  Of 
course,  the  disclaimer  must  be  added  that  the  calculations  only  indicate  the 
overall  effect,  and  that  it  is  certainly  possible  that  some  error  in  k^  has  been 
offset  by  a corresponding  error  in  k^. 

This  assignment  of  k^  was  based  exclusively  upon  comparisons  of  the  exponen- 

18 

tial  decay  of  CH2O  in  the  mixtures  with  concentrations  near  4.5  x 10  molecules/ 

3 

cm  . Additional  calculations,  using  the  rate  constants  of  C3,  were  done  at  2.3  x 
18  3 

10  molecule/cm  . Values  of  k.  obtained  for  the  1%  mixtures  here  were  in  good 

O 

agreement  with  those  at  the  higher  concentrations;  this  was  consistent  with  the 
experimental  observations.  A final  set  of  calculations  was  used  to  check  the 
entire  [CH20]-time  profile.  Here  the  purpose  was  to  use  the  C3  rate  constants, 
which  gave  good  fits  to  k^,  to  see  if  these  values  were  also  able  to  describe 

a 

the  early  time  behavior  in  the  experiments.  At  the  higher  temperatures  and 
correspondingly  short  reaction  times,  the  observed  emission  signal  is  affected 
by  the  finite  slit  width  as  well  as  the  response  time  of  the  detector.  To  ac- 
count for  these  factors  the  calculated  concentration-time  profile  was  first 
corrected  for  the  1 ys  detector  response  time;  this  modified  concentration-time 
profile  was  then  integrated  over  a trapezoidal  slit  function.  The  detector  efficiency 


!! 

(j 

i Mas  assumed  to  be  unity  for  the  1.5  itm  slit  width  and  then  drop  linearly  to  zero 

ij  at  a distance  2.0  mm  beyond  each  silt  edge.  Here  to  was  defined  as  the  time 

at  which  the  reflected  shock  was  first  visible  to  the  detector.  Typically,  to 

i 

was  5 ps  earlier  than  the  arrival  of  the  shock  front  at  the  silt  midpoint.  Re- 
suits  of  such  a calculation  are  shown  In  Fig.  1.  The  fit  appears  to  be  quite 
reasonable;  note  that  a shift  of  the  calculated  values  by  lys  would  noticeably 
Improve  the  fit.  Such  a shift  Is  compatible  with  the  estimated  1-2  ys  uncer- 
tainty In  locating  to.  Fits  comparable  to  that  shown  In  Fig.  1 were  also  ob- 
served at  the  higher  temperatures.  At  lower  temperatures,  the  reaction  times 
^ were  sufficiently  long  that  silt  effects  can  be  neglected.  It  Is  In  this  region 

that  the  data  suggest  an  Induction  period  prior  to  the  onset  of  exponential  de- 
cay. Comparison  of  profiles  In  this  region  again  Indicated  good  agreement.  Par- 
ticularly satisfying  was  the  fact  that  the  calculations  showed  the  same  variation 
of  Induction  time  with  mixture  dilution  that  was  seen  experimentally.  In  summary, 
the  rate  constant  combination  C3  appears  to  accurately  model  all  aspects  of  the 
observed  CHgO  decay.  In  retrospect.  It  would  appear  that  the  failure  of  Eq.  [A] 
to  account  for  various  aspects  of  the  data  was  not  due  to  a deficiency  In  the 
mechanism;  the  problem  was  obviously  related  to  the  mechanistic  analysis  via 
the  steady-state  assumption. 

A comparison  of  the  rate  constant  for  formaldehyde  dissociation  to  the  total 
rate  constant  for  decay  gives  one  an  estimate  of  the  Importance  of  the  secondary 
reactions  In  formaldehyde  pyrolysis.  The  ratio  k^/k^  ranges  from  1%  for  1% 
mixtures  at  1800  K to  'v*  40%  for  0.1%  mixtures  at  2500  K.  This  variation  is  also 
mirrored  In  the  product  distribution.  At  the  lower  temperatures,  the  calcula- 
tions Indicated  that  [Hgl/CH]  » 1 once  there  had  been  any  appreciable  decay  of 
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CH2O.  At  2500  K this  ratio  is  somewhat  higher  than  unity  for  1%  mixtures,  but 
becomes  much  less  than  unity  for  the  0.1%  mixture,  it  is  interesting  to  note 
that  the  product  distribution  from  the  chain  decomposition  mechanism  may  often 
resemble  that  expected  if  the  dissociation  were  to  go  directly  to  molecular  pro- 

I • 

ducts,  i.e.,  reaction  (1  ). 

The  observed  rates  of  formaldehyde  decay  in  all  of  the  shock  tube  experi- 
ments suggest  that  the  ' rate  constant  assignment  from  the  CH^/02  flame  work^ 

is  much  too  large.  At  1800  K,  the  shock  tube  results  indicate  a total  decay 

153  1-1  12 

rate  constant  of  2 x 10"  cm  molecule"  s'  , whereas  k^ ' * 2 x 10"  at  this 

temperature.  Furthermore,  the  shock  tube  work  suggests  that  secondary  reactions 

I 

play  a major  role  in  the  decay  scheme.  Thus,  decay  of  CH2O  via  reaction  (1  ) 
must  have  a rate  constant  much  lower  than  k.,  and  the  present  work  suggests  that 

o 

I 

the  total  decay  can  be  adequately  modeled  without  inclusion  of  reaction  (1  ). 

It  would  appear  that  the  molecular  dissociation  path  is  not  a significant  reac- 
tion channel  in  the  high  temperature  pyrolysis  of  formaldehyde. 
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TABLE  I 


SUMMARY  OF  EXPERIMENTAL  OBSERVATIONS 


(a) 

Mixture 

A 

1.01%  CH^O 


B 

1.01%  CH2O 

C 

1.00%  CH^O 


D 

0.50%  CH^O 


L 


T/K 

18  -3^^^ 

M/10  **  molecule  cm  ^ 

1^/10”^®  cm^  molecule 

2545 

4.77 

1010 

2350 

4.62 

426 

2175 

4.48 

201 

2065 

4.37 

121 

1970 

4.28 

77.3 

1945 

4.24 

49.8 

1885 

4.21 

• 32.8 

1795 

4.11 

20.0 

2710 

2.45 

1630 

2700 

2.45 

1180 

2585 

2.42 

1620 

2540 

2.40 

792 

2470 

2.40 

812 

2385 

2.33 

700 

2275 

2.28 

391 

2220 

2.31 

302 

2510 

4.73 

791 

2465 

4.69 

595 

2330 

4.59 

410 

2250 

4.54 

359 

2150 

4.46 

156 

747 

564 

568 

385 

271 

332 

178 

115 

46.4 
32.7 
21.0 

22.5 
8.77 


2010  4.26  54.5 

2005  4.24  49.7 

2005  4.27  48.8 


2505 

4.74 

2445 

4.59 

2390 

4.68 

2385 

4.68 

2280 

4.62 

2270 

4.51 

2220 

4.50 

2055 

4.33 

1970 

4.33 

1925 

4.19 

1830 

4.16 

1810 

4.17 

1700 

3.98 

A 


1 

TABLE  I . (Continued) 

(a) 

Mixture  T/K 

18  -3^^^ 

M/10  ° molecule  cm  ^ 

_g  3 _i  ) 

k /lO  ° cm"*  molecule  s 

E 2435 

4.62 

— — 

740 

2380 

4.58 

434 

0.51X  CH9O  2255 

4.46 

197 

^ 2160 

4.41 

151 

2040 

4.25 

. 76.5 

1960 

4.19 

42.0 

1845 

4.08 

26.5 

1765 

4.01 

12.2 

0.1  ox  CH2O 


2595 

4.63 

410 

2480 

4.57 

444 

2475 

4.54 

379 

2350 

4.47 

257 

2300 

4.45 

188 

2220 

4.37 

87.9 

2190 

4.30 

73.7 

2085 

4.30 

39.3 

2030 

4.19 

28.2 

1930 

4.11 

17.1 

1855 

3.95 

10.8 

1790 

3.95 

5.57 

^*^Balance  of  all  mixtures  was  argon. 
^**Hotal  concentration. 

^^^ka  = -(d  Jin  (Signal  )/dt)/M. 


U(A.B,C)  -20.69  0.267  197.4  4.7  0.16 
.0.5%  (D.E)  -20.24  0.601  208.3  10.2  0.19 
0.1%  (F)  -20.18  0.495  222.1  6.8  0.20 


TABLE  III.  COMPARISON  OF  CALCULATED  AND  OBSERVED 

KAi  iirr  Ar  l,  (1) 


VALUES  OF  k. 


0.1%  CH«0 

0.5%  CHaO 

1. 

0%  CH„0 

1800K 

2500K 

1800K 

2500K 

1800K 

" 2500K 

(2) 

>erved'  ' 

6.09  ± 1.20 

390  ± 77 

14.7  ± 2.8 

725  ± 136. 

19.4  ± 3 

.1  778  ± 124 

Cl 

26.9 

676 

46.6 

1140 

51.0 

1440 

C2 

6.57 

548 

6.60 

594 

6.49 

578 

C3 

7.07 

349 

14.6 

652 

16.2 

808 

C4 

6.27 

294 

13.9 

594 

15.2 

752 

RATE  COHSTANTS  FOR  CALCULAT I Or.'S 


k,t3) 


Set 

In  A 

la 

In  A 

la 

In  A 

ia 

In  A 

Cl 

-16.3 

301 

-23.9 

16.0 

-22.1 

61.4 

.-21.8 

C2 

-16.3 

301 

-23.9 

16.0 

-25.5 

79.4 

-21.8 

C3 

-14.3 

364 

-23.9 

16.0 

-22.1 

61.4 

-21.8 

C4 

-16.1 

339 

-22.4 

37.6 

-22.1 

61.4 

-21.8 

(1) 

Units 

-1  fi  3 

are  10"  cm  molecule 

s"\ 

• 

(2) 

Error 

bars  are  ±a. 

4 

(3) 

3 

A in  cm  molecule' 

s-1  E 

5 » ‘'a 

In  kJ/mol . 

LEGENDS 


Fig.  1 A typical  experimental  plot  of  log  (IR  signal)  vs.  time.  To  is  defined 
as  the  time  when  the  reflected  shock  first  enters  the  field  of  view  of 
the  detector.  Mixture  B,  T * 2150  K,  M = 4.46  x'lO^®  molecule/cm^. 

O = observed  points  (for  clarity  every  other  data  point  has  been 
omitted);  • = calculated  profile  (see  text).  Arrows  indicate  re- 

gion used  to  determine  decay  rate;  solid  line  is  a least-squares  fit 
to  the  observed  data. 

Fig.  2 Arrhenius  plot  of  for  those  shocks  with  M 'v  4.5  x 10^®  molecule/cm^. 

O = mixtures  A,  B,  and  C.  (1%  CH20/Ar);  • = Mixtures  D and 

E (0.5%  CH20/Ar);  ^ = Mixture  F (0.1%  CH20/Ar).  Lines  are  least- 

square  fits;  error  bars  shown  are  one  standard  deviation. 
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These  data  were  analyzed  In  terms  of  the  mechanism; 

CH2O  + M = HCO  + H + M (1] 

CH2O  + H = HCO  + H2  (2] 

HCO  + M = H + CO  + M (3; 

HCO  + H = H2  + CO  (4] 

Using  literature  values  for  k2»  ks,  and  k4,  a reasonable  fit  to  a11 
the  data,  both  In  terms  of  temperature  and  concentration  dependence,  was 
obtained  with  ki  = 6.0  x 10”’exp(-364  kJ/RT)  cm®  molecule"®  s"®.  Evidence 
Is  also  presented  to  suggest  that  the  reaction 

CH2O  + M = CO  + H2  + M (r 

does  not  play  a significant  role  In  the  high  temperature  pyrolysis  of 
formaldehyde. 
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